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Abstract 

Background: Estrogen receptors play an important role in mediating estrogen action on target tissues, and the 
estrogen is relevant to male infertility. Single nucleotide polymorphisms (SNPs) in estrogen receptors may be 
associated with the risk of male infertility. A variety of case control studies have been published evaluating this 
association. However, the accumulated studies have shown inconsistent conclusions. 

Methods: To further determine the potential association between the four common SNPs (rs2234693, rs9340799, 
rsl 256049 and rs4986938) in estrogen receptors gene and male infertility, this meta-analysis was performed 
according to the 10 published case control studies. The odds ratio (OR) and 95% confidence interval (CI) were used 
to evaluate the strength of the associations. 

Results: It was revealed that the sub-group analysis by the ethnicity, for the rs2234693, a significant association in 
the comparison of CC vs. TT (OR = 0.61, 95% CI: 0.40-0.93), CT vs. TT (OR = 0.67, 95% CI: 0.49-0.93) and CC + CT vs. TT 
(OR = 0.66, 95% CI: 0.49-0.89) in the Asian population with male infertility. For rs9340799 polymorphism, increased 
risks were observed for the comparison of AA vs. GG (OR = 1 .75, 95% CI: 1 .1 5-2.68) and AA vs. GA + GG (OR = 1 .38, 
95% CI: 1.02-1.88). For rsl 256049 polymorphism, the comparison of the GA vs. GG (OR = 1.52, 95% CI: 1.00-2.31) 
and AA -i- GA vs. GG (OR = 1 .74, 95% CI: 1 .03-2.94), also increased risks present in Asian and Caucasian population, 
respectively. 

Conclusions: The rs2234693C allele was associated with the decreased risk for male infertility; however, the 
rs9340799AA genotype and the rsl256049GA genotype were associated with an increased risk for male infertility. 

Keywords: Male infertility. Polymorphisms, Estrogen receptors 



Background 

Male infertility is an important cause of couple's inability 
to bear children in 20% to 25% of total cases and the eti- 
ology of nearly half of the cases remains idiopathic [1,2]. 
Approximately 15% of male infertile cases, genetic factors, 
including chromosomal aberrations and single gene mu- 
tations, may result in spermatogenic failure and sperm 
dysfunction [3,4]. The traditional view of estradiol as the 
'female' hormone and of testosterone as the 'male' hor- 
mone has been challenged due to the increased interest in 
elucidating the role of estrogen in males [5]. Estrogens are 
produced in the male reproductive system by Sertoli cells, 
Leydig cells, and germ cells [6,7]. In addition, studies re- 
vealed that estrogens reduce testosterone production from 
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Leydig cells and reduce Sertoli cell numbers in adult when 
they are given during development [8,9]. The estrogens 
can also disrupt fetal Leydig cell development, inhibit 
apoptosis of human postmeiotic germ cells, and increase 
spermatogonial number per testis [8-12]. The physical 
functions of estrogens were involved in the estrogen re- 
ceptors (ERs). Moreover, ERs are members of the nuclear 
receptor (NR) superfamUy that mediates the pleiotropic 
effects of estrogen in a diverse range of developmental 
and physiological processes, playing an important role in 
mediating estrogen action on target tissues [13,14]. 

ERs have been identified to be two subtypes of ERa 
and ERp. ERa is a 595-amino acid protein [15] encoded 
by the ERsl gene on chromosome 6q25, and ERp is a 
530-amino acid protein [16] which encoded by the ERs2 
gene on chromosome 14q22-24 [17]. Genetic screening 
for the ERa gene locus has revealed several polymorphic 
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sites [18], and two polymorphisms located in ERa intron 
1(T/C transition, rs2234693) and in 50 bp downstream of 
the former one (G/A transition, rs9340799) have been 
widely concerned. In addition, the ER|3 genes have been 
described with two silent G/A polymorphisms (rsl256049 
and rs4986938) [19]. To date, epidemiological studies have 
been carried out to evaluate the association between ER 
polymorphisms and male infertility. However, the results 
remain inconsistent (Table 1) [5,7,19-26]. In order to get a 
more precise estimation of the association between poly- 
morphisms in ERs and risk of male infertility, this meta- 
analysis was performed based on ten eligible previously 
published studies. 

Methods 

Identification and eligibility of studies 

To identify all articles that examined the association of 
ERs polymorphisms with male infertility, a comprehen- 
sive systematic bibliographic search through the medical 
databases PUBMED, attempting to cover all medical 
papers published between 1950 and 2013, using the fol- 
lowing keywords and subject terms: "male infertility", 
"polymorphism" and "estrogen receptors" or "ERs". The 
synonyms of polymorphism (rs2234693, rs9340799, 
rsl256049, and rs4986938) were also used as the key- 
words in the search. The studies were excluded if they 
were not English language papers or human subject. 
References in retrieved articles were screened in which 
case reports, meta-analyses and review articles were ex- 
cluded. In addition, studies were identified by a manual 
search of the references lists of reviews and retrieved 
studies. All the studies were included if they met the 
following criteria: (I) about the rs2234693, rs9340799, 
rsl256049, and rs4986938 polymorphisms and male 
infertility, (II) from a case control study, (III) geno- 
type frequencies could be derived. The reasons for 
exclusion of articles were listed in the Additional file 1: 
Text SI. 

Data extraction 

Two authors (Tian-Fu Li and Qiu-Yue Wu) extracted all 
data independently that met the inclusion criteria and 
reached the consensus for any controversy. The main 
characteristics of the enrolled studies were listed in the 
Table 1, including: (I) the first author's last name, (II) 
year of publication, (III) ethnicity, (IV) source of control 
groups (population- or hospital-based controls), (V) 
genotyping methods, (VI) the polymorphism sites, (VII) 
characteristics of studies, (VIII) Case/Control counts, 
(IX) Hardy- Weinberg equilibrium in the controls. Data 
were extracted separately for each ethnic groups catego- 
rized as Caucasian and Asian. However, no African was 
identified in this study. 



Statistical analysis 

The risk of male infertility associated with the four poly- 
morphisms of the ERs gene was estimated for each study 
by odds ratio (OR), together with its 95% confidence 
interval (CI), respectively. The four polymorphisms were 
evaluated for the associations with male infertility sus- 
ceptibility based on four genetic models. To contrast, 
the wild-type homozygote (WW), we first estimated the 
risk of the rare allele homozygote (RR) and heterozygous 
(WR) genotypes on infertility, then evaluated the risk of 
infertility under a dominant model (RR + WR vs. WW). 
In addition, recessive model associations were also esti- 
mated (RR vs. WR + WW). Moreover, stratified analyses 
were also performed by ethnicity (Asian and Caucasian). 
The statistical significance of the pooled OR was deter- 
mined with the Z-test and a P-value of <0.05 was con- 
sidered significant. Heterogeneity across the studies was 
evaluated by Chi-square test based on Q test [27] and 
was considered significant if P <0.05.A fixed-effect mo- 
del using the Mantel-Haenszel method and a random- 
effects model using the DerSimonian and Laird method 
were used to pool the results [28]. In addition, the fixed- 
effect model was used as well when there was no hetero- 
geneity across results of the studies, or the random-effect 
model. Moreover, a sensitivity analysis, by which a single 
study in the meta-analysis was deleted each time to deter- 
mine the influence of the individual data set to the overall 
pooled OR, was performed to assess the stability of the re- 
sults. To test the publication bias. Funnel plots and Egger's 
linear regression test were applied [29]. Hardy- Weinberg 
equilibrium in the controls of each study was calculated 
using a web-based program [30]. All statistical tests for 
this meta-analysis were performed with STATA version 
10.0 (Stata Corporation College Station, TX, USA). 

Results 

Characteristics of studies 

A total of 10 eligible case control studies with the publi- 
cation dates ranged from 2002 to 2013 met the prespeci- 
fied inclusion criteria (shown in the Figure 1), including 
five studies of Asian population [19-21,23,24] and five 
studies of Caucasian population [5,7,22,25,26]. To deter- 
mine the SNPs, two different genotyping methods such 
as PCR-RFLP [5,19-26] and TaqMan assays [7] were ap- 
plied. All subjects were received comprehensive andro- 
logical examination, and the patients were divided into 
three types: oligozoospermia (sperm count <20 x lO^/mV), 
azoospermia and oligoasthenoteratozoospermic (OAT). 
The studies' exclusion criteria and inclusion criteria were 
listed in the Table 1. In addition, the sources of controls in 
these studies were mainly population-based. The distri- 
bution of genotypes in the controls of all studies was con- 
sistent with Hardy- Weinberg equilibrium except for the 



Table 1 Summary of published studies included 



Author Year Race Source Method Polymorphism Characteristics of study patients Case/control counts HWE 

of control sites (Control) 



Meng [19] 


2013 


Asian 


PB 


PCR-RFLP 


rs2234693, 
rs9340799, 
rs 1256049, 
rs4986938 


Age: 25-38 years (mean age 32.1 ±5.2 years). Exclusion criteria: 
abnormal karyotypes, deletions of the Y chiromosome, orchitis, 

1 j_ 1 ■ 1 ■ -J. 1 1 ■ 1 j_ 1 1 £ j-\ 

varicocele, cryptorchidism, congenital bilateral absence or the vas 
deferens, hypogonadotropic hypogonadism, and iatrogenic 
infertility. 


"17:83/82, CT:96/126, CC:25/44; AA:151/ 
1 48, AG:42/89,GG:1 1/15; GG: 1 03/1 27, 
AG:91/102, AA:10/23; GG:155/193, 
AG:41/48, AA:8/1 1 


0.712, 0.793, 
0.699, 0.001 


Zaiata [5] 


2013 


Caucasian 


PB 


PCR-RFLP 


rs2234693, 
rs9340799 


Inclusion criteria: same ethnic origin (Caucasians). Exclusion 
criteria: varicocele, hormonal therapy, hypogonadism, smoking, Y 
chromosome deletions and karyotype abnormalities. The ages of 
were not shown in the article. 


TT:33/14, CT:32/27, CC:16/19; AA:28/8, 
AG:36/32, GG: 1 7/20 


0.468, 0.389 


Ogata [20] 


2012 


Asian 


PB 


PCR-RFLP 


rs 1256049 


Age: 32-52 years (median 41.0 years). Inclusion criteria: no 
extragenital anomalies, seminal tract obstruction, varicocele, Y 
chromosomal microdeletion, or retrograde ejaculation; normal 
karyotypes. 


GG:68/64, AG:49/45, AA:8/1 0 


0.604 


Bianco [7] 


2011 


Caucasian 


PB 


TaqMan 
assays 


rs2234693, 
rs9340799, 
rs 1256049 
rs4986938 


Age: 36.1 ±6.5 years. Exclusion criteria: chromosome anomalies, 
azoospermia factor (AZF) microdeletions, smoking, alcoholism, 
occupation, varicocele, and cryptorchidism. 


TT:30/37, CT93/1 1 1, CC:64/68; AA:80/ 
100, AG:79/88, GG:20/28; GG:1 72/201, 
AG:15/15, AA:0/0; GG:43/28, AG:60/103, 
AA:84/85 


0.468, 0.221, 
0.597, 0.712 


Safarinejad [21] 


2010 


Asian 


PB 


PCR-RFLP 


rs2234693, 
rs9340799, 
rs 1256049, 
rs4986938 


Age: 31. 6 ±4.8 years (range 25-40 years). Inclusion criterion: two 
years with no reason for their infertility. Exclusion criteria: 
varicocele or testicular torsion, urinary tract infections, 
endocrinopathy, karyotype anomalies, Y-chromosome 
microdeletions, use of drugs, leukocytospermia, a BMI of 
30 kg/m2 or greater. 


Tr:49/33, CT70/86, CC:45/45; AA:62/41, 
AG:77/95, GG:25/28; GG:142/152, 
AG:21/8, AA:1/4; GG:65/80, AG:82/63, 
A A: 17/21 


0.486, 0.034, 
0.000, 0.132 


Lazaros [22] 


2010 


Caucasian 


PB 


PCR-RFLP 


rs2234693, 
rs9340799, 
rs 1256049, 
rs4986938 


Age: 33.2 ±67.5 years. Exclusion criteria: hypogonadotropic 
hypogonadism, obstructive syndromes of the seminal tract, 
microdeletions of the Y chromosome, karyotypic abnormalities. 


rr:6/20, CT:14/40, CC:9/25; AA:5/13, 
AG:13/43, GG:1 1/29; GG:26/80, AG:3/5, 
AA:0/0; GG:7/1 7, AG:12/36, AA:10/32 


0.609, 0.652, 
0.779, 0.246 


Khattri [23] 


2007 


Asian 


PB 


PCR-RFLP 


rs 1256049 


Age: 23.24 ±2.06 years. Exclusion criteria: obstruction, 
endocrinological defect, injuries, karyotypic abnormality, 
Y-chromosome microdeletions. 


GG:397/231, AG:46/21, AA:0/0 


0.490 


Omrani [24] 


2005 


Asian 


PB 


PCR-RFLP 


rs 1256049 
rs4986938 


Exclusion criteria: genetic causes of infertility, such as Klinefelter 
syndrome or Ychromosome microdeletions. The ages of patients 
were no shown in the article. 


GG:103/194, AG:17/9, AA:0/1 ; GG:51/86, 
AG:57/88,AA:12/30 


0.023, 0.339 


Aschim [25] 


2005 


Caucasian 


PB 


PCR-RFLP 


rs 1256049, 
rs4986938 


Exclusion criteria: Klinefelter syndrome or Y-chromosome 
microdeletions, a history of cryptorchidism were excluded. 
The ages of patients were no shown in the article. 


GG:92/177, AG:14/8, AA:0/1 ; GG:1 1/82, 
AG:48/79, AA:47/25 


0.015, 0.394 


Kukuvitis [26] 


2002 


Caucasian 


PB 


PCR-RFLP 


rs2234693, 
rs9340799 


Exclusion criteria: any known aetiologies (varicocele, infections of 

accessory glands, cryptorchidism, homozygous form of 
(3-thalassemia). The ages of patients were no shown in the article. 


TT:38/18, CT:38/25, CC:33/21 ; AA:30/1 0, 
AG:45/28, GG:34/26 


0.083, 0.594 



PB, Population Based; PCR-RFLP, Polymerase Chain Reaction-restriction Fragment Length Polymorphism; HWE, Hardy-Weinberg equilibrium; BMI, body mass index. 
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24 of records identified through PubMed searching 



6 records removed after screening the abstract (not ER 
poljTnorphism and male infertilit}', rev!e\\') 



18 of full- text articles assessed for eligibility' 



7 of full-text articles excluded (2 not case-control studies, 
4 other SNP of ER, 1 article without independent data) 



11 of full-text articles assessed for eligibility^ 
articles assessed 



1 article excluded (without available data) 



10 anicles included in analysis 



Figure 1 Flow chart of studies identified with inclusion and exclusion criteria. 



study [19] in rs4986938, the study [21] in rs9340799 
and rsl256049, study [24] in rs 1256049 and study 
[25] in rsl256049, which were tested in the sensitivity 
analyses. 

Quantitative synthesis 

Wide variation of four polymorphisms allele frequencies 
across different ethnicities was observed. For rs2234693, 
the frequency of T allele was 53.13% (95% CI: 49.74- 
56.52) in the Asian controls, which was higher than that 
in Caucasian controls 44.82% (95% CI: 41.48-48.16) as 
shown in Figure 2A. For rs9340799, the frequency of G 
allele in the Asian controls (32.45%, 95% CI: 29.27-35.63) 
was lower than that in Caucasian controls (46.71%, 
95% CI: 43.36-50.06) as shown in Figure 2B. In Figure 2C, 
we could find that the frequency of G allele for the 
rsl256049 in the Asian controls (87.34%, 95% CI: 
85.88-88.81) was lower than which in Caucasian con- 
trols (96.92%, 95% CI: 95.82-98.02). In contrast, the 
frequency of G allele in Caucasian controls (48.46%, 
95% CI: 45.32-51.60) was lower than that in Asian group 
(73.39%, 95% CI: 70.93-75.85) for the rs4986938 in 
Figure 2D. 

For rs2234693 polymorphism, significant differences 
were observed for the comparison of CC vs. TT, CT vs. 
TT and CC -i- CT vs. TT. Sub-group analysis by the eth- 
nicity revealed a significant association in the com- 
parison of CC vs. TT (OR = 0.61, 95% CI: 0.40-0.93, 



^'heterogeneity = 0.670), CT VS. TT (OR = 0.67, 95% CI: 
0.49-0.93, Pheterogeneity = 0.358), CC + CT VS. TT (OR = 

0.66, 95% CI: 0.49-0.89, Pheterogeneity = 0.593) and C alleles 

vs. T alleles (OR = 0.78, 95% CI: 0.64-0.96, ^heterogeneity = 

0.681) in the Asian population, as summarized in 
Table 2. 

For rs9340799 polymorphism, increased risks were 
observed for the comparison of AA vs. GG and AA 
vs. GA -I- GG. Sub-group analysis by ethnicity revealed 
increased risks (AA vs. GG: OR = 1.75, 95% CI: 1.15- 

2.68, ^heterogeneity = 0.174; AA VS. GA + GG: OR = 1.38, 

95% CI: 1.02-1.88, ^heterogeneity = 0.062) in the Caucasian 
population, also for the AA vs. GA -i- GG and A alleles vs. 
G alleles, a significant association was observed in Asian 

population (OR = 1.93, 95% CI: 1.42-2.62, ^heterogeneity = 
0.768; OR =1.49, 95% CI: 1.18-1.87, ^heterogeneity = 0.375) 

as summarized in Table 3. 

For rsl256049 polymorphism, significant differences 
were observed for the comparison of GA vs. GG, AA -i- 
GA vs. GG and AA vs. GA -i- GG. For the comparison of 
the GA vs. GG, AA -i- GA vs. GG, increased risks present in 
Asian and Caucasian population, respectively (GA vs. GG: 
OR = 1.52, 95% CI: 1.00-2.31, ^heterogeneity = 0.038; AA-i- 

GA vs. GG: OR = 1.74, 95% CI: 1.03-2.94, Pheterogeneity = 

0.275). All data were concluded in the Table 4. In contrast, 
a decreased risk was also observed for the comparison AA 

vs. GA-hGG (OR = 0.55, 95% CI: 0.31-0.97, Pheterogeneity = 

0.818) in Asian population. For the rs4986938, there was 
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Figure 2 The allel frequencies of the four polymorphisms in the controls may vary by ethnicity. (A) rs2234693, (B) rs9340799, (C) rs 1256049, 
(D) rs4986938. Star or dot denotes outliers. 



no significant association observed in all comparisons (data 
were not shown). 

Test of heterogeneity 

Among the four polymorphisms, a significant hetero- 
geneity was apparent among heterozygote comparison 

for the rsl256049 (GA vs. GG: ^heterogeneity = 0.047) 

(Figure 3). Two studies [7,19] were identified to con- 
tribute to substantial heterogeneity, and it was decreased 
when the study was removed respectively {P = 0.065, 
P = 0.075). Sensitivity analysis revealed that the two inde- 
pendent studies [7,23] were the main cause of heterogen- 
eity for the rsl256049. Heterogeneity was decreased when 
these studies were removed (GA vs. GG: ^heterogeneity = 
0.320, f = 14.7%). Although the genotype distributions in 
four studies did not follow Hardy- Weinberg equilibrium, 
the corresponding pooled ORs were not materially altered 
by excluding the studies. 



Publication bias 

Begg's funnel plot and Egger's test were performed to as- 
sess the publication bias of the currently available litera- 
ture. For the rs2234693, rs9340799 and rs4986938, the 
shape of the funnel plots did not reveal any evidence of 
obvious asymmetry in all comparison models. Then, the 
Egger's test was used to provide statistical evidence for 
funnel plot symmetry. The results also did not show any 
evidence of publication bias. However, for the rsl256049, 
as shown in the Figure 4, the shape of the funnel plots 
seemed asymmetrical in the heterozygote and dominant 
comparisons, suggesting the presence of publication bias. 
Then, the Egger's tests were adopted to provide statis- 
tical evidence of funnel plot asymmetry. As expected, 
the results showed obvious evidence of publication bias 
{t = 2.53, P = 0.044 for GA vs. GG; t = 2.71, P = 0.035 for 
AA + GA vs. GG). To adjust for this bias, a trim-and- 
fill method developed by Duval and Tweedie [31] was 



Table 2 Stratification analyses of genetic susceptibility of rs2234693 polymorphism to male infertility 



Category 


Cases/ 
controls 


CC vs. 


TT 




CT vs. 1 


rr 




CC-l-CT vs. TT 




CC vs. CT + TT 




C allele vs. T allele 




OR(95% CI) 


pa 




OR(95% CI) 






OR(95% CI) 






OR(95% CI) 


pa 




OR(95% CI) 






Total 


774/841 


0.72(0.54-0.96) 


0.282 


20.1 


0.74(0.58-0.94) 


0.533 


0 


0.73(0.58-0.91) 


0.368 


7.5 


0.90(0.71-1.13) 


0.465 


0 


0.84(0.71-1.01) 


0.190 


32.8 


RACE 


































Asian 


368/41 6 


0.61(0.40-0.93) 


0.670 


0 


0.67(0.49-0.93) 


0.358 


0 


0.66(0.49-0.89) 


0.593 


0 


0.83(0.58-1.18) 


0.257 


22.1 


0.78(0.64-0.96) 


0.681 


0 


Caucasian 


406/425 


0.83(046-1 .23) 


0.175 


39.5 


0.83(0.57-1.19) 


0.460 


0 


0.81(0.58-1.14) 


0.232 


30 


0.95(0.70-1.29) 


0.395 


0 


0.87(0.64-1.18) 


0.103 


51.5 



0-25, no heterogeneity; 25-50, modest, 
value of Q-test for heterogeneity test. 
Bold numbers mean statistically significant results. 



Table 3 Stratification analyses of genetic susceptibility of rs9340799 polymorphism to male infertility 



Category 


Cases/ 
controls 


AA vs. GG 




GA vs. GG 




AA -1- GA vs. GG 




AA vs. GA -1- GG 




A allele vs. G allele 




0R(9S% CI) 






OR(95% CI) 






OR(95% CI) 


pa 


F 


OR(95% CI) 






OR(95% CI) 






Total 


774/841 


1.67(1.21-2.32) 


0.392 


3.9 


1.03(0.76-1.39) 


0.764 


0 


1.27(0.96-1.68) 


0.796 


0 


1.63(1.32-2.03) 


0.077 


49.7 


1.39(1.13-1.68) 


0.172 


35.3 


RACE 


































Asian 


358/416 


1 .56(0.93-2.62) 


0.714 


0 


0.81(0.49-1.34) 


0.524 


0 


1.13(0.71-1.82) 


0.952 


0 


1.93(1.42-2.62) 


0.768 


0 


1.49(1.18-1.86) 


0.375 


0 


Caucasian 


406/425 


1.75(1.15-2.68) 


0.174 


39.7 


1.17(0.81-1.71) 


0.847 


0 


1.35(0.95-1.92) 


0.569 


0 


1.38(1.02-1.88) 


0.062 


59.2 


1.39(0.97-1.81) 


0.109 


50.5 



value of Q-test for heterogeneity test. 

0-25, no heterogeneity; 25-50, modest heterogeneity; 50, high heterogeneity. 
Bold numbers mean statistically significant results. 



Table 4 Stratification analyses of genetic susceptibility of rs1 2S6049 polymorphism to male infertility 



Category 


Cases/Controls 


AA vs. GG 




GA vs. GG 




AA -1- GA vs. GG 




AA vs. GA -1- GG 




A allele vs. G allele 


OR(95% CI) 






OR(95% CI) 






OR(95% CI) 






OR(95% CI) P= 




OR(95% CI) 


pa ,2 


Total 


1378/1478 


0.57(0.32-1.01) 


0.940 


0 


1.59(1.12-2.25)'' 


0.047 


50.9 


1.30(1.05-1.61) 


0.075 


45.7 


0.55(0.32-0.96) 0.920 


0 


1.29(0.97-1.72) 


0.068 46.9 


RACE 






























Asian 


1056/991 


0.57(0.32-1.01) 


0.853 


0 


1.52(1.00-2.31)'' 


0.038 


60.6 


1.23(0.98-1.56) 


0.064 


55 


0.55(0.31-0.97) 0.818 


0 


1.19(0.86-1.65) 


0.067 54.4 


Caucasian 


322/487 


0.64(0.03-1 5.86) 






1 .87(0.92-3.80) 


0.208 


36.2 


1.74(1.03-2.94) 


0.275 


22.6 


0.58(0.02-1438) 




1 .66(0.99-2.77) 


0.372 0 



value of Q-test for heterogeneity test. 
'^Random-effects model was used when a P value, 0.05 for heterogeneity test; otherwise, fixed-effects model was used. 

0-25, no heterogeneity; 25-50, modest heterogeneity; 50, high heterogeneity. 
Bold numbers mean statistically significant results. 



Li et al. Reproductive Biology and Endocrinology 2014, 12;79 
http://www.rbej.eom/content/1 2/1 /79 



Page 9 of 13 



Stujy 
ID 



OR (95% CI) 



Weight 



Meng 2013 
Ogata 2012 
Bianco 2011 
Safarineisd 2010 
Lazaros2010 
Khattri2007 
Omrani200E 
Aschim200E 

Oi'erall (l-squared = SO.9%, p = 0.047} 

NOTE: Weights ars Ifom rantjom effects analysis 
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1.10(0.75, 1.61) 
1.02(0.60, 1.74) 
1.17(0.56,2.46) 
2.81 (1.21,6.55) 
1.85(0.41,8.26) 
1.27(0.74.2.19) 
3.56(1.53,8.26) 
3.37(1.36,832) 
1.59(1.12,2.25) 



20.24 

16.56 

12.06 

10.38 

4.49 

16.29 

10.43 

9 54 

100.00 



8-32 



Figure 3 Forest plot for the overall association between rsl 256049 polymorphism and male infertility for random effects. For GA vs. GG 

each study was shown by the point estimate of the OR (the size of the square is proportional to the weight of each study) and 95% Ci for the OR 
{extending lines); the pooled OR and 95% CI were shown by diamonds. 



implemented. Trimming was based on fixed-effects mo- 
del, and the adjusted estimates obtained by using the ran- 
dom effects model were ORs of 1.17 (0.78-1.74) for 
GA vs. GG and 1.08 (0.75-1.54) for AA+GA vs. GG 
in the Figure 5. Although Meta-analysis with or with- 
out the trim-and-fill method also ends up with same 
conclusions, but the ORs were not statistically signifi- 
cant difference. So it was indicated that the results of 
these studies were not statistically robust. 

Discussion 

The present meta-analysis, including 1568 cases and 1602 
controls from 10 case control studies, explored the associ- 
ation between the ERs polymorphisms and male infertility. 
The results indicated that rs2234693C allele was associ- 
ated with decreased risk of the male infertility, particularly 
in the Asian population. In contrast, rs9340799AA geno- 
type was observed as a risk factor for infertility in both 
Asian and Caucasian population, and rsl256049GA ge- 
notype was associated with an increased risk for develop- 
ing male sterility. However, the rs4986938 polymorphism 
was not associated with male infertility. In addition, we 
tried to find the data in the available database, such 



as PUBMED [32], National Human Genome Research 
Institute GWAS Catalog [33] and GWAS Central [34], but 
we found no relevant genome-wide association (GWAS) 
study about these four polymorphisms. 

Estradiol has been reported as a survival factor for 
germ cells [11], involving in the induction of oxidative 
DNA damage, and the aberrant level of estrogen may 
lead to impaired sperm production [35-37]. It has been 
shown that free radicals inhibit steriodogenesis by inter- 
fering with cholesterol transport to the mitochondria 
and/or the catalytic function of P450 enzymes, which 
leads to an increase in lipid per oxidation and decline in 
the antioxidant barrier [38]. Moreover, estrogens can regu- 
late mitochondrial function by increasing nuclear respira- 
tory factor-1 (NRF-1) expression [5]. Specifically, estradiol 
stimulates mitochondrial function through a genomic me- 
chanism of ER action involving direct ERa and ER|3 inter- 
action with an oestrogen response element in the NRF-1 
promoter [39]. In vivo knockdown experiments have in- 
dicated that estradiol stimulates NRF-1 transcription and 
consequently increases mitochondrial biogenesis through 
ERa activity but not through ER|3 activity in MCF-7 breast 
cells [40] . This findings indicates that ERa polymorphisms 
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can increase mitochondrial activity via NRF-1 transcrip- 
tion in human ejaculated spermatozoa, presenting them 
with high motility [22]. 

The mechanisms behind altered ERs function in sub- 
jects with polymorphisms remain unclear. The polymor- 
phism rsl256049 located at the splice acceptor site just 
prior to exon 8 in ERp [41] and may potentially affect 
the splicing of this exon, leading to proteins with differ- 
ent properties than the wild-type ERp [42,43]. In addi- 
tion, studies have reported the polymorphism could also 
have a direct effect through changing the nucleotide se- 
quence and thereby the secondary structure of the ER(3 
mRNA, possibly leading to changes the function of mRNA 
[44,45]. It has been reported that ERa gene polymor- 
phisms (rs2234693 and rs9340799) may modalate the ef- 
fect of oestradiol on CYP19, which encodes aromatase 
expression, disrupting the gene causes a decline in sperm 
numbers and loss of male infertility [46,47]. 



The precise role of estrogen receptors in male fertility 
status is understood. Some findings suggest that specific 
polymorphisms of the ERa, and ER|3 genes which confer 
a lower sex hormone binding globulin (SHBG) and thus 
a stronger unbound estrogen effect, may adversely affect 
human spermatogenesis [48,49]. SHBG is involved in 
both delivering reproductive hormones to target tissues 
and controlling the concentration of androgens and es- 
trogens in the serum and tissues [50]. Pavlovich et al. 
[51] demonstrated that infertile men with severe oligo- 
zoospermia had significantly lower T (testosterone) and 
higher E2 (estradiol) concentrations than fertile control 
subjects, resulting in an elevated T/E2 ratio. 

Identifying the source of heterogeneity is one of the 
most important goals of the meta-analysis. Thus, we stra- 
tified the studies only according to ethnicity (because the 
sources of the controls were selected through population- 
based, and the method used was the only one different). 
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adjusted estimates obtained by using the random effects model for GA vs. GG and AA + GA vs. GG. 



Stratified analysis by ethnicity revealed that there was no 
difference between the European population and Asian 
population, suggesting that different ethnicities and en- 
vironmental exposures may have no influence on the sus- 
ceptibility of male infertility, and more studies should be 
accumulated to reveal the difference. In addition, for the 
rsl256049, sensitivity analysis revealed that the three inde- 
pendent studies [7,22,23] were the main source of hetero- 
geneity. Heterogeneity was decreased when these studies 
were removed. For these three studies, the sample size 
was not sufficient and the numbers of rsl256049AA geno- 
type was both zero. These two points may be the main 
reason for the heterogeneity in the performed analysis. 
For the rsl256049, there was obvious evidence of publica- 
tion bias. As the same with heterogeneity, the numbers of 
the cases and controls of the wild-type homozygote in 
these three studies [7,22,23] were too small to keep the 



results statistically robust, so it maybe the key factor for 
the bias. Using a proper and representative subject is very 
important in reducing bias in such genotype association 
studies. 

There are still some limitations in this meta-analysis. 
Firstly, there were only ten literatures enrolled in this 
meta-analysis, the sample size was not big enough to have 
substantial power exploring the real association. Secondly, 
the detailed information (such as life-style, age, and work) 
could not be traced, so that our unadjusted estimates 
should be confirmed by further studies. In addition, an in- 
dividual with a clinical disorder was not a result of the sin- 
gle gene that is disrupted, but that the genetic disruption 
was embedded within the context of that individual's en- 
tire genome and environment exposure [52] . In fact, some 
other genes related to fertility could also play an important 
role in spermatogenesis. 
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Conclusions 

In summary, this meta-analysis suggested that the 
rs2234693C allele was the protective factor for male 
infertility, the rs9340799AA genotype was associated with 
an increased risk for infertility, and the rsl256049GA 
genotype was also the negative factor. 
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